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ABSTRACT 

We use Hubble Space Telescope (HST) near-infrared imaging to explore the shapes of the surface brightness 
profiles of bulges of SO-Sbc galaxies at high resolution. Modeling extends to the outer bulge via bulge-disk 
decompositions of combined HST - ground based profiles. Compact, central unresolved components similar to 
those reported by others are found in ^84% of the sample. We also detect a moderate frequency (^34%) of 
nuclear components with exponential profiles which may be disks or bars. Adopting the Sersic r'/" functional 
form for the bulge, none of the bulges have an r'^"* behaviour; derived Sersic shape-indices are < n >= L7 ±0.7. 
For the same sample, fits to NIR ground-based profiles yield Sersic indices up to n = 4-6. The high-n of ground- 
based profiles are a result of nuclear point sources blending with the bulge extended light due to seeing. The low 
Sersic indices are not expected from merger violent relaxation, and argue against significant merger growth for 
most bulges. 



Subject headings: galaxies:spiral — galaxies: structure - 
galaxies:photometry 

1. INTRODUCTION 

Whether galaxy mergers (e.g. Kauffmann et al. 1996) or in- 
stead disk instabilities (Pfenniger & Norman 1990) may be 
viewed as the dominant process of formation and growth of the 
central bulges in spiral and SO galaxies depends to some extent 
on the shape of the surface brightness profiles of bulges. Simu- 
lations of mergers of similar-size galaxies yield spheroidal sys- 
tems following the r'/'* law (e.g. Barnes 1988). The accretion 
of dense satellites onto bulge-disk galaxies also leads to r'/'* 
bulges (Aguen-i, Balcells & Peletier 200L herafter ABPOl). In 
contrast, disk instabilities lead to bulges with approximately ex- 
ponential profiles (Combes et al. 1990). 

Over the past decade several studies have stressed that bulges 
of late- to intermediate-type spiral galaxies are best fit with pro- 
files of exponential to r''^ shapes (Andredakis & Sanders 1994; 
de Jong 1996; Courteau, de Jong & Broeils 1996; Seigar & 
James 1998; Carollo, StiavelH & Mach 1998; MacArthur et al. 
2002). However, many bulge-disk decompositions are still per- 
formed using r'/^ bulges only, especially at high-z, e.g. Simard 
et al. (2002); the r'/"* law is often cited for massive or early-type 
bulges: Carollo and collaborators (e.g. Seigar et al. 2002 and 
references therein) routinely classify bulges into r'/"* and expo- 
nential types. Do massive bulges follow the r'/"^ model? Do 
bulges come in two flavors, perhaps pointing at massive bulges 
growing by early mergers and less-massive bulges growing by 
disk instabilities? 

Fitting bulge profiles with the Sersic model fj,(r) ~ r'/" pro- 
vides an avenue for measuring the true shapes and concentra- 
tions of the light distribution of bulges. Various studies that in- 
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elude early-type bulges (Andredakis, Peletier & Balcells 1995, 
hereafter APB95; Khosroshahi et al. 2000; Graham 2001 ; M61- 
lenhoff & Heidt 2001) find a continuous distribution of Sersic 
shape indices n that scales with bulge luminosity from n < 1 to 
n > 4. At the high-« end, these results depend strongly on the 
innermost regions of the profiles. Inner disks, compact sources 
and star formation, present in galaxy nuclei (e.g. Phillips et al. 
1996; Carollo et al. 2002; Rest et al. 2001; Pizzela et al. 2002), 
contribute light that cannot be disentangled from the spheroid 
light at ground-based resolutions (Ravindranath et al. 2001). 
Extra central light results in higher-n profile fits. 

Nineteen galaxies from the Balcells & Peletier (1994, BP94) 
sample were imaged with HST/NICMOS at F160W (Peletier 
et al. 1999, Paper I). We have performed bulge-disk decom- 
positions, and analyze profiles, global parameter relations and 
nuclear cusp slopes in a companion paper (Balcells et al. 2002, 
Paper III). In this Letter, we focus on the shape indices obtained 
from Sersic fits to the bulge profiles. The high resolution of 
the HST data allows one to disentangle nuclear disks/bars and 
compact sources from the extended spheroid light. Modeling 
these components, we study Sersic indices for the bulge pro- 
files and compare them to those obtained from ground-based 
profiles. We characterize the magnitude distribution of point 
sources found in the galaxy centers and discuss the relevance 
of low-n profiles for bulge formation models. Preliminary re- 
sults on nuclear cusps are presented in Balcells et al. (2001). 

A key ingredient of the analysis of bulge profiles, overlooked 
in the HST-hased studies of bulges cited above, is the model- 
ing of the extended bulge and disk light distributions when fit- 
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ting analytic functions to the bulge light. Here, composite pro- 
files linking HST profiles with ground-based K-b&nd profiles 
are used so that bulge-disk decompositions can be performed. 
We assume a Hubble constant of Hd = 75 km s"^ Mpc"^ . 

2. DATA 

We work with a subset of the BP94 diameter-limited sample 
of inclined, early-to-intermediate type disk galaxies classified 
as unbarred in the UGC (Nilson 1973); see Paper I and refer- 
ences therein. None has a Seyfert or starburst nucleus. Ground- 
based surface brightness profiles in the /T-band for the BP94 
sample, derived from ellipse fits to UKIRT/IRCAM3 images, 
are given in Peletier & Balcells (1997), while two-dimensional 
bulge-disk decomposition and Sersic fits to the /T-band bulge 
profiles are published in APB95. 

The subsample studied here comprises 19 galaxies of types 
SO to Sbc that were imaged with NICMOS on HST (camera 2, 
19"xl9", 0.075 arcsec/pixel). Data reduction is described in 
Paper I, while details of the data analysis are given in Paper 
III. We derive elliptically-averaged surface brightness profiles 
and isophotal shapes, from 0.075"to typically 10", using the 
galphot package (J0rgensen et al. 1992). We keep the cen- 
ters fixed and allow the elUpticity and position angles of the 
isophotes to vary. Linking the HST-F160W and ground-based 
K profiles is done following Paper I, with details given in Paper 

m. 

3. BULGE-DISK FITS 

The combined HST plus ground-based profiles, corrected 
for foreground Galactic extinction (Schlegel et al. 1998), 
(1 + z)'^ cosmological dimming and A'-correction, were fit- 
ted with a PSF-convolved Sersic -plus-exponential model using 
the code from Graham (2001) modified to use a Moffat PSF 
(FWHM=0.19"[NICMOS Handbook], /?=2.5). No internal ex- 
tinction corrections were applied. A quasi-Newton algorithm 
(Kahaner, Moler & Nash 1989) was used to converge on the 
solution resulting in the smallest rm.s. scatter, with 5 free pa- 
rameters (disk /Lto, h; bulge iJ,e,Re,n). Each fit was individually 
inspected. When an obvious local minimum was found, the 
code was run again with new initial conditions; when the ap- 
parent global minimum was found, the code was also again run 
with different initial conditions to provide confirmation. The 
final fits for every galaxy are shown in Paper III. 

Such Sersic-exponential bulge-disk fits, covering the entire 
radial range down to 0.075", tend to show strong residuals, up 
to 0.3 mag, over the entire bulge. A typical example is shown in 
Figure la. The residual profile shows a characteristic wave pat- 
tern. This is caused by a strong central Ught excess, as inferred 
from fits excluding the central 1-2 arcsec. We then run the fit- 
ting code including either an additional central delta function 
point source (PS); a central exponential; a point source plus an 
exponential; or a central Gaussian component with freeFWHM 
to the fitting function. In all cases each component is PSF- 
convolved before fitting. For each galaxy we select the simplest 
fit, i.e. minimum number of components, for which the residual 
profile falls below an amplitude of 0. 1 mag/(arcsec)^. Figure lb 
shows the profile from Figure la fitted with an additional cen- 
tral PS, which yields an excellent fit. The total light in the inner 
1-2" can also be approximated with power laws (Balcells et al. 
2001; Boker et al. 2002). The Nuker model also fits the in- 
ner few arcsec of the early-type galaxy profiles well (Byun et 
al. 1996). However, these approaches break down further out 



(e.g. Lauer et al. 1995), do not allow a bulge-disk decomposi- 
tion, and are not discussed further here. Our fits extend over 2.5 
decades in radius, and describe both bulge and disk. 
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Fig. 1. — (a) The combined /f5r+ground-based H-bmd surface 
brightness profile of NGC 5443 (Sb), with PSF-convolved Sersic bulge 
and exponential disk components (soUd lines). The residuals (data mi- 
nus model) are shown in the lower panel, (b) The same surface bright- 
ness profile fitted with a PSF-convolved S6rsic law, an exponential 
disk, plus a central point source (dashed line). Residuals are shown in 
the lower panel. 

Of the 19 galaxies, 12 yield a satisfactory fit with the addition 
of a single PS, 2 with the addition of a single nuclear exponen- 
tial, 4 with a PS plus an exponential, and one with a pure Sersic 
bulge. The frequency of PS is thus 84 ± 8%, and that of in- 
ner exponentials is 32± 11%. Only one galaxy (NGC 5577, 
Sbc) does not require an unresolved source and can be fit with 
a simple Sersic bulge plus exponential disk; however, the pro- 
file for this galaxy has lower S/N than the others owing to its 
lower surface brightness. Most galaxies with PSs yield a good 
fit with a Gaussian inner component instead; in those cases, the 
FWHM of the inner Gaussian, before PSF convolution, is gen- 
erally well below the PSF FWHM, recovering the PS solution. 
Some galaxies with inner exponential components also admit a 
good fit with a marginally resolved Gaussian instead. Because 
these galaxies show high ellipticity and pointy isophotes in the 
nuclear region, we choose the PS-plus-exponential fit instead of 
the Gaussian inner component. Nuclear isophotes are described 
in Paper HI. 

The Sersic index for the bulges ranges from 0.5 to 3.0, 
with a mean of <n>= 1.7 ±0.7. The range is significantly 
lower than that obtained by APB95 using ground-based data 
alone, which reaches n = 6: Figure 2b shows that is sys- 
tematically lower than righ, especially for galaxies with ngb>3. 
None of our bulges reaches the de Vaucouleurs r'/'* behaviour. 
It appears that fits to ground-based profiles reach Sersic indices 
n > 4 because the Ught from i/5r-unresolved central sources, 
plus in some instances nuclear disks or bars, when convolved 
with typical ground-based seeing, link smoothly with the the 
extended bulge profile and mimic higher-n Sersic profiles. 

Figures 2c,d plot the distribution of rihst vs. the bulge-to- 
disk luminosity ratio (B/D) and bulge A'-band absolute mag- 
nitude, respectively. B/D are derived from the best-fit param- 
eters. Bulge absolute magnitudes are derived from the galaxy 
/T-band magnitude given in APB95 and our computed B/D. As 
noted by APB95, some of the galaxy magnitudes are proba- 
bly slightly underestimated when the galaxy overfills the frame. 
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Fig. 2. — Dependence of the i/ST-resolved bulge Sersic index rii,,, on global properties of the parent galaxy: (a) vs galaxy morphological type index T. (b) 
ni,s, vs the bulge Sersic index Ugi, obtained from ground-based data (APB95). (c) iihst against the bulge-to-disk ratio, (d) Bulge index n^, against the bulge K-bicad 
absolute magnitude. Filled circles: bulges, this work. Crosses: bulges from the de Jong & van der Kruit (1994) sample, as analyzed by Graham (2001) (cf. his 
Figure 14). 



Figures 2c,d show that log(n/„() correlates with log(B/D) and 
with the bulge absolute A'-band magnitude (Spearman rank- 
order correlation coefficients -0.52 and -0.51, significances 
97.8% and 97.4%, respectively). Thus, the correlation of Sersic 
index n with B/D and bulge luminosity found in ground-based 
data (APB95; Graham 2001) remains valid when fits including 
central PS components are used. We compare our HST-hased 
n distributions to those obtained by Graham (2001) using the 
sample of spirals of de Jong (1996) (Figure 2c,d, cross sym- 
bols). The loci occupied by the two samples are similar. A 
trend with morphological type T (Figure 2a) is apparent when 
our distribution is extended to later types by adding values from 
the Graham (2001) analysis. 

4. POINT SOURCE PROPERTIES 

We showed in §3 that, for 16 out of the 19 galaxies stud- 
ied, adding a central PS improves the fit to the observed pro- 
file. Absolute //-band magnitudes for the PS are in the range 
-14 > Mh.ps > -18. The distribution of Mhj>s with morpho- 
logical type T is quite flat (Figure 3a). Mh,ps correlates with 
log(n) and with Mk.buI (Spearman rank-order correlation co- 
efficients of -0.63 and 0.77, significances 99.6% and 99.96%, 
respectively), see Figures 3b,c. The line in Figure 3c is an or- 
thogonal regression to the Mh,ps—Mk,buI relation, which yields 

Lk,ps/Lk,q = \0''\LKfiui/lO'^LK,Qf-'^^''\ (1) 

using a constant color H-K = 0.23 for the nuclei (Persson, Fro- 
gel & Aaronson 1979). 

Nuclear sources have been reported in previous studies of 
early-type galaxy profiles based on HST imaging. CaroUo et al. 
(2002) detect nuclear sources in 43 of their 69 late-type galax- 
ies, with absolute magnitudes -10 > Mh.ps > -17.7. The bright 
end is similar to ours. The fainter luminosities are consistent 
with the expected low bulge luminosity of their later type galax- 
ies, although the lack of a bulge-disk decomposition prevents 
us from comparing theirs and our trends of PS luminosity with 
spheroid luminosity. Ravindranath et al. (2001), using two- 
dimensional fits to HST/NICMOS F160W images of an E/SO 
sample, report unresolved central components for ~50% of 
their targets. Figure 3d shows their distribution of point source 
magnitude vs. absolute B-band galaxy magnitude together with 
our distribution of PS magnitude vs. fi-band absolute bulge 
magnitude; the latter was computed using the Bt galaxy mag- 
nitudes from the RC3 and the bulge-disk ratio derived from the 



//5r+ground-based profile fits. Their point source magnitudes 
are sUghtly fainter than ours. Their fainter PS magnitudes and 
their lower detection rate may be consequences of their choice 
of the Nuker model vs. our use of the Sersic model to fit the un- 
derlying bulge light distribution; the Nuker law is steeper than 
the Sersic law near the origin, and diverges at r = 0, whereas the 
Sersic law has a finite value of /i(0). Also, using galaxy magni- 
tudes for their sample instead of bulge magnitudes shifts their 
points to the right in the diagram. 

5. DISCUSSION 

Measuring the correct value of the Sersic index n for bulges 
is important for several reasons. Overestimating n leads to 
an overestimate of the bulge and luminosity. Trujillo et al. 
(2001) show that fitting an r^/^ profile with a r'/"^ function 
yields an overestimate of rg by a factor of 3-5 and of the bulge 
luminosity by a factor of about 2. The luminosity ratio B/D 
is further overestimated as more of the light in the bulge-disk 
transition region is assigned to the bulge. The converse occurs 
when n is underestimated, eg. when exponential fits are forced 
onto n > I Sersic profiles (Graham 2001). Thus, small lumi- 
nosity components (typically L^j>s/La:.b«/ ~ 10"^ **, eqn. 1) may 
significantly alter the derived global parameters of bulges. Test 
fits excluding inner profile regions suggest that, in the absence 
of //>Sr imaging, the effects of central components on the bulge 
shape index n are greatly diminished by excluding the central 
FWHM of the profile, as done by Stiavelli et al. (2001) for 
dwarf ellipticals. Point sources should present less of a prob- 
lem for late-type spirals given the low Sersic indices, low bulge 
luminosities and hence low PS luminosities expected from Fig- 
ure 3c. 

The nuclear PSs, amounting to a fraction of order lO"-' of 
the bulge luminosities, are easy to accomodate within either a 
merger origin or a secular evolution origin for bulges. The ob- 
served scaling of PS and bulge luminosities need not indicate 
an internal origin for the PS as nuclear gas deposition during 
a merger may scale with the masses of the merging objects. If 
bulges are evolved bars, the problem of growing a nuclear star 
cluster is similar to that of feeding an active galactic nucleus 
(Maciejewski et al. 2002). 

Given the persistent description of bulges as r'/'* spheroids 
in galaxy formation models, it is worth pointing out that none 
of flie bulges foUow flie r^/^ law. APB95, de Jong (1996), 
Graham (2001) and Mac Arthur et al. (2002) have shown that 
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Fig. 3. — Dependence of central point source //-band absolute magnitude ps on global properties of the parent galaxy: (a) pg vs. galaxy morphological 
type index T. (b) Mu ps vs. Sersic index n of the bulge, (c) M/j pg vs. //-band absolute bulge magnitude. The line is an orthogonal regression to the distribution. 
(d) Mh,ps vs. absolute /?-band magnitude. Filled circles: our bulges; absolute magnitude is that of the bulge. Crosses: SO galaxies from Ravindranath et al. (2001). 
Open circles: E galaxies from Ravindranath et al. (2001). Total galaxy absolute magnitudes are used for SO and E galaxies. 



late- to intermediate-type bulges have r^/" profiles in the range 
0-5 ^ n ^ 2.0. Our work establishes a similar range, 0.5 < 
n < 3.0, for intermediate- to early-type bulges i.e. those which 
are still generally described with the r'/"* model. Extrapolating 
from the rising envelope of the n - Mk relation in Figure 2d, 
only bulges with < -26 may reach n = 4. The r^/"* law 
is interpreted as the consequence of violent relaxation during 
a clumpy collapse (van Albada 1982), and is obtained after 
mergers of similar-size disk galaxies (Barnes 1988). Simula- 
tions of collisionless satellite accretion onto bulge-disk galax- 
ies further show that bulges evolve smoothly from exponential 
profiles to r'/^ profiles when doubUng their mass due to accre- 
tion events (ABPOl). At face value then, our distribution of 
Sersic shape indices n (Figure 2c,d) suggests that the vast ma- 
jority of bulges did not acquire their present structural prop- 
erties by growing through collisionless mergers. Kinematic 
data might show whether, besides a low Sersic index n, these 
bulges show other properties typical of disks (Kormendy 1993; 



Falcon-Barroso et al. 2002). 

Have deviations from the r^^^ law in A^-body models gone 
undetected? Sersic fitting is only starting to be applied in A^- 
body studies. Hjorth & Madsen (1995) show that violent relax- 
ation in a finite volume, combined with subsequent escape of 
positive-energy stars, leads to a distribution function with de- 
viations from the r'/"* law; from their figures, a shallow central 
potential might possibly fit n < 4 Sersic profiles. Accretion of 
low-mass satellites might heat the bulge only moderately and 
deposite the compact satellite nucleus at the center of the rem- 
nant. This could lead to a n ~ 1-2 bulge with a central bright- 
ness spike not unlike what is observed. But very long merger 
times and slow growth rate for low-mass satellites are a prob- 
lem for this hypothesis. Recent cosmological SPH simulations 
of galaxy formation (Scannapieco & Tissera 2002) find n Ri 1 
for bulges resulting from secular evolution, while merging rises 
n to about n = 4, again pointing at low-n profiles being related 
to non-merger formation histories. 
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